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Abstract. A Jordan- Wigner (JW) transformation for the t — J and 
Hubbard models is described. Introduced are holon and doublon par- 
ticles for hole and double occupied sites. There is only a single spin 
sector particle. Flux tubes occur in a naturel fashion, within a specific 
gauge, when the method is adapted to two dimensions. In order to ac- 
commodate three dimensions "flux sheets" are defined. The adaptation 
of the method to the t — J-model in the context of high Tc is described. 

While the Jordan- Wigner transformation [1] is well known in the context spin only 
models such as the Heisenberg oi X — Y model it does not seem to be widely known 
that a similar approach can be used for doped systems such as the t — J or Hubbard 
models and that this provides a useful alternative formulation of the two dimensional 
such models in the context of high Tc superconductors. The purpose of the present 
Letter is to describe such a formalism along with some elementary conclusions which 
are evident in the basis of this formalism. 

Recall first the basic JW transformation in one dimension. It is trivial for 5=1/2 
that obey on-site Fermi commutations rules and that, 

fl^e^-^ll\'^-at (1) 

obeys such rules even when the site indices are different. Here Qn = fnfn^ i-^-i ^^1^ 
number of particles at site n. It is then straightforward to show [1] that, e.g., the one 
dimensional Heisenberg Hamiltonian Ti. = J'^n ' '^n+i becomes, 

+ H.C.) + J ^ fUnfi+lL+l - ^^T (2) 
n n 

where, with iV| — J2n fnfm the last term reflects an effective chemical potential fj, — J. 

This can, in fact, be generahsed to include charge. Since cr+ is considered to be a 
creation operator it is implicit that the vacuum is the down spin ferromagnetic state. 
Creating a /^-fermion implies an up spin at the site n. In the usual way a 6jj-boson 
corresponds to a true vacuum state |0) without any particle, i.e., to the presence of a 
hole, while djj-boson implies that the site n is in the state cj^jcj^lO), which is doubly 
occupied. The physical operators are given by. 



a+=e--i::".'3'"/t, (3a) 
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c\n = {flK + die-'- e^-*^ , (36) 

cl = (^'ne- - e"^T , (3c) 

where now Qn = .fnfn^Hi^n^did^^ counts all particles and where there is a constraint 
that Qn < 1, i.e., the particles have a hard core. It is straightforward to check that 
these obey the appropriate commutation rules when the site indices are different. 

When the the site indices are the same it is not the case that {c^^^, c|^} = 1 if the 
algebra of the auxiliary particles is applied without strictly applying the constraint. 
What is easily verified is that the matrix elements of say on the physical sub-space 

are correctly given, however it is implied that in the product, e.g., c|„C|„ it is necessary 
to include a physical complete set of states between the two operators. In this regard 
the present scheme differs from the traditional auxiliary particle scheme [2] for which 
this is not necessary. It is implied that care must be exercised when constructing, in 
particular, the auxiliary particle version of the Hamiltonian. Every physical operator, 
for the site n, can be written in terms of the product of one or two of the new auxiliary 
particles since any operator can be written in terms of: 

U)ai = i-(4/n+&l&n+4rf„), iT)(Ti = 4/„ io)(oi = 6t6„, \n){n\ = did^ 
|o)(Ti = 6U„. in)(Ti = 4/„, in)(o| = 4^. iT)Ui = e-^'S:-«-/t, 

|0)(i I = e-'-^l=\ ^-e-^Tfot , I ||)(| I = e-^-Sr=\ ^-e^^^rft, (4) 

and their Hermitian conjugates. 

For the one dimensional f — J- model the phase operators, e"^^T and e 

m=l 

cancel and the result is, 

W = J2(fnKbUjn+l + + H.C.) + i Y.(fnfn+1 + H.C.) 

n n 

+ J E(24/n + bibJ{2fl^Jn+i + bUiK+i) - (5) 

n 

The Hubbard model becomes: 

n = -t^{flbJn+ifn+i + dlfji+idn+, + 6t6„+i + dUn+i + H.C.) + Udld^ 

n 

n 

+4/n&Uie-^'^^-^'5^ +6t4+i/„+ie-^'E::\'3'^ +F.C.). (6) 

Even in one dimension, for the Hubbard model, the JW formulation would seem 

not very useful because of the presence of the factors of e^™ 5^m=i . The situation is 
even worse in two dimensions. The transformation can be used for a two dimensional 
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Fig. 1 (a) The zig-zag path used for the Jordan- Wigner transformation, (b) The 
"string" flux lines used to evaluate the sign changes due to the Jordan- Wigner trans- 
formation. 



system once a mapping to one dimension is defined. A possible such map, shown 
in figure la, is to use a zig-zag path whence even for the t — J-model the factors of 

e-"i:::\<3"> appear. Each site is ordered by a single site label and the sum in this 
factor simply counts the number of particles with smaller such labels. It is well known 
that, e.g., the up spins of the problem, which are equivalent of hard core bosons, 
might be converted to fermions by attaching a unite flux tube. The awkward factors 

of e"*'^^™=i "^"^ in fact represent such flux tubes a certain "string" gauge. 

Consider the appropriate such phase factor for vertical exchange J-coupling, be- 
tween sites n and m, as shown in figure lb. The appropriate term in the Hamiltonian 
is easily seen to be, 



•ErJ+i«V(/i/„ + i/.c.); n<m, 



(7) 



where the sum X^^n+i Qm counts the number of particle between the two sites foUow- 

ing the prescribed path. The resulting phase factor e ^q=,^+l " can be evaluated 
using the diagramme method illustrated in this figure. To each particle are attached 
two lines or "strings" , as shown, and which leave the system without crossing the 
zig-zag path, i.e., leave by the open ends. It follows that, 



(8) 



where L is the number of lines which cross a straight line joining the sites n and n' . Put 
differently, L is the number of lines cut when the particle hops from site n to site m. In 
a fermion plus flux tube formulation the same phase factor would be exp[— i a - dr\, 
where a is an appropriately defined vector potential. The identification is therefore: 



a ■ dr — ttL. 



Any path which encloses a single particle cuts two lines and has. 



a ■ dr = 2TT, 



(9) 



(10) 
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corresponding to a single flux quantum. Evideirtly different mappings from one to 
two dimensions corresponds to a different version of this JW string gauge. The phase 
factors appropriate to correlation function are also determined, e.g., 

([5+(t),5-(0)]) = (-l)^-([4(i),/„(0)]), (11) 

where L„„i is the number of lines crossed for any path joining the sites n and m which 
does not pass through any particle. (It need not necessarily follow bonds.) 

With this method, e.g., the full t — J- model for a two dimensional square lattice 

is, 

^ = - E (tfih^jfj + tiMbj + H.c.) + ^ ^(///,. + H.c.) 

<ij> <ij> 

+ i E + blhm^jfj + b]bj) - J(27Vt + No), (12) 

<ij> 

where tij = {—1)^'H and Jij = {—l)^'^J and where Lij is the number of strings 
crossed in going for the site i to j as described above. Here Nq = J2i ^J^i- There is 

no factor (— l)^*-) in the term since the factors e *'^^m=n+i cancel. (The 

term involves the exchange of particles and hence each crosses the same number of 
strings and the total number is necessarily even.) These factors also cancel in static 
interaction term with the prefactor j. (Here there is no movement of particles and 
hence no strings are crossed.) 

It is rather obvious, for the, a; = 0, undoped case, that for the physical correla- 
tion function ([5+(i),S-(0)]) = (-l)^'""([/t(t), /„(0)]) the string prefactor (-1)^"- 
cancels against the similar factors in Jij contained in the propagator and that hence 
this quantity is gauge invariant. In the doped case the same result is less obvious 
since the spin particle can propagate by an exchange with a holon (or doublon) and 
such a hopping process does not "see" the flux. However, the holon (or doublon) 
must move on a closed path and in the end the (—1)^"™ factors involved in closing 
the path compensate. Similarly for charge propagators, e.g., within the t — J-model 
{'^■rcl^{T)c^^{0)) = (— l)'^'""(rT-6„(T)6j„(0)) and similar statements can be made. 
Thus in general, e.g., 

([5+(i),5-(0)])=expH / a.rfr1([/t(t),/„(0)]). (13) 

J n 

When making potential gauge transformations it is worth noting that the sense of a 
unite flux tube is immaterial. In particular, for a bipartite lattice, alternating the 
sense of the flux according to the sub-lattice is an easy way to avoid making gauge 
transformations which create (unphysical) currents at the boundaries. 

Important also is the observation that, with the JW gauges, all matrix elements 
of Ti. are real. For such a Hamiltonian matrix it is a trivial fact that either the ground 
state vector is (i) real (or can be made real with a simple change of phase) and therefore 
carries not currents or (ii) is degenerate. In case (ii) the real and imaginary parts of 
the ground state vector correspond to degenerate states which also carry no current. 
Case (ii) does not exclude the possibility of broken symmetry ground states which 
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Fig. 2 (a) The 7r/4 string gauge transformation. There is an effective half fliix quantum 
in the shaded squares, (b) The equivalent with 7r/8 strings. 

carry currents. If the degeneracy of the ground state is associated with the spin part 
of the wave function then there must be a finite value of S{S +1) ~ 5^, i.e., such 
degenerate ground states must have a net ferromagnetic moment, albeit as small as 
S=l/2. 

Up to this point the flux tubes have been attached to the particles. It is equally 
possible to attach the flux tubes, and the false charge which sees them, to the "empty" , 
i.e., down spin sites. This scheme will be assumed in what follows. 

It is clearly possible to associate an operator with a flux tube at a given site. In 
the operator u^{{Qi}) — exp[— i ^ a-dr\ the integral from the origin to site n passes 
by any path which does not include particles and specifically along some path which 
does not pass through any sites. The argument {Qi} indicates that u„ is a function of 
the position of all of the particles, or more precisely reflects the position of all of the 
down spin sites since it is these which determine a. Clearly this is a unitary operator 
such that w]^u„ — w„u,]j — 1. For a JW gauge itj^ — it„. Trivially these ujj operators 
commute with each other but anti-commute with the particle operators. This formal 
development permits the, e.g., Ujblbj term to be written as: 

-tblu.ulb^. (14) 

and e.g., — u^b.^. It should be noted that these unitary operators always appear in 

such a fashion that m| is associated with the creation of an down spin site along with 
a flux tube while only occurs when such a tube is destroyed. Using these operators 
the t — J model, 

^ = E (f^b^b]f, + bhhu]b^ + H.C.) + ^{.f^,U,ulf^ + H.C.) 

<ij> <ij> 

+ 1 E (2///. + b\bm]f, + blb^) - + No). (15) 

<ij> 

The present JW-gauge does not reflect the symmetries of the lattice and in ad- 
dition the interaction between fermions which is induced by the tt-Aux strings is of 
infinite range. The usual gauge with a = {l/r)(t) is an evident alternative but is 
not particularly adapted to lattice problems. It is simpler to adapt a different string 
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W sheet 



Fig. 3 (a) A 18 site bi-layer is viewed as being to a single folded plane and the zig-zag 
path 1,2,..., 17, 18 is designated on that plane. (Site 10 is hidden behind the sheet 
but can be seen in (b) and (c).) The JW-sheet is atatchcd to site 11 which lies on the 
path between the illustrated bond between sites 5 and 14. That the sheet cuts this 
bond implies that the sign of the bond changes, (b) A similar JW-sheet attached to 
site 13 also lies on the path between these sites and the sheet cuts the bond idicating 
again that such a particle causes a change of sign. However in (c) the sheet is attached 
to site 15 is not on the path between 5 and 14 and indeed the sheet does not cut the 
bond indicating that there is no change in sign. In general there are many sheets as in 
(d). The sheet attcahed to site 4 in the lower plane does not cut the bond and there is 
therefore only a single change of sign associated with the particle at site 11. Explicity 
the sheet which passes throught site 11, coordinates x = —1, y = 0, and z = 1, was 
generated by z = 1 - 0.4tan"i[3(y + 4tan~i{3(a; + 1)}]. 

gauge. Each tt string can be thought of as made of four ±7r/4 strings, where the 
sign is arbitrary. Separating one such string and moving it so it cuts different bonds, 
see figure 2a, amounts to some gauge transformation. (The flux pattern is such that 
there is an effective half flux tube in each of the two shaded plaquettcs.) If instead 
7r/8 strings are used it is possible to have the gauge convention shown in figure 2b. 
Continuing this bifurcation process will lead to interactions which ~ (1/r). The arrow 
indicates the sign of the phase change. As in electrostatics, a positive magnetic flux 
tube with arrows heading away from the tube will have positive sign changes when a 
particle crosses a line heading in a clockwise direction. 

It is of interest to generalise the present approach to bi-layers and three dimensions 
in general. While the fermion-plus-flux-tube type of aynons do not generalise to three 
dimensions there is no real problem with applying the JW transformation. A bi-layer 
can be considered as a single sheet which is folded over and the single sheet is then 
reduced to one dimensions using again, e.g., the zig-zag path of figure la. Fully three 
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dimensional system can similarly be; mapped to a corrngatcd sheet and then the sheet 
reduced to one dimension. On two dimensions the JW-transformation is effected by 
attaching tt flux strings to the particles. The generalisation to three dimensions is 
to attach sheets to each particle. Consider, by way of illustration the 18-site bi-layer 
shown in figure 3. By design, if a sheet is cut at the level of a given plane then the 
resulting line has the same topology as the flux strings of flgure 1 and the sign of a 
intra-plane bond (not shown) is given, as before, by (—1)^ where L is the number of 
such lines, or now sheets, which cross a bond. Each site on the zig-zag path, indicated 
by the black lines, has be labled and advances in the sense 1, 2, . . . , 17, 18 as shown in 
figure 3a. Focus attention on the white bond which connects sites 5 and 14. Within 
the JW transfromation the sign associated with this bond is (—1)'^ where S is the 
number of particles on the sites 6, 7, ... , 12, 13 which lie between the ends of the bond 
on the prescribed path. The number S is counted using the JW-sheets. In figure 3a 
the sheet is attached to a particle at site 11, identified by an arrow. This sheet cuts the 
bond indicating, correctly, that this site must be included in S. Similarly, figure 3b, 
the sheet attached to site 13 cuts the bond and is to be reflected in S. However site 
15 which lies on the same row as this latter site should not have an effect on S and 
indeed, figure 3c, the sheet attached to this site does not cut the bond. Of course, in 
general there are many particles and many sheets which potentially are reflected in 
the factor (—1)"^. In figure 3d are shown two sheets, the one attached to site 4 is not 
in the sequence 6, 7, ... , 12, 13 and the sheet does not cut the bond between 5 and 14, 
while, as described above that attached to the site at 11 does, and should, cut the 
bond. Notice the sheets never cross the zig-zag path. 

When two particles are interchanged it is necessarily the case that one particle 
hops through the JW-sheet of the other particle an odd number of times thereby 
converting the fermions to hard core bosons. Once defined with a specific zig-zag 
path, deforming the sheets amounts to a gauge transformation, and a tt sheet can 
be separated into n, Ti/n sheets which can be deformed almost at will. In this way 
the interactions between particles, reflected by the sheets, can be made short ranged. 
(Also, clearly, substituting a ^-sheet for the 7r-sheet permits an extrapolation between 
fermions and hard core bosons in three dimensions.) 

Whatever the dimension, the hare vacuum \)-n/2 reflects the absence of par- 
ticles and comprises the spin state with Sz = —N/2. It is important that the 
bare ferromagnetic vacuum state is highly degenerate. A new vacuum |)_Ar/2+i = 

(5'(^/M_jJy2 )l)-JV/2+i is obtained by acting with the total spin raising operator 

where = y/S{S -|- 1) — nm. The physical spin vacuum IS") comprises the Fermi 
spin sea. To be specific define this as, 

\s) = n fi\)-N/2, (16) 

where eg oc -7^; 7g = cos kx -\- cos ky. Given that there are /-particles (and that 
this is an integer) this is a spin singlet. Although this description of IS') is appealing, it 
is sometimes more appropriate to imagine \S) as derived from the Sz = bare vacuum 
Dq . In the wave function fl is replaced by p^ = \/2(| + Snz) and the empty (down) 

sites are projected out using — V2{^ — Snz)- With this \S) = f ({5„^})|)o, where 
{Snz} is the collection of local z-component spin operators. 
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The states f^\S) and / -|S') are momentum k excitations with Sz — ±1 respec- 
tively and are related by a particle-hole symmetry at half-filling. In addition to these 
/-particle excitations, are similar states generated by the flux tube operators u^. A 
arbitrary state with a given vacuum Sz = nis of the form Fi{{Snz})\)n ^^'^ typical 
matrix element of is 

= J\F^u,ifjf, + fjl)F,\)„ 

n{\F2fiFl\)n+l^m,n + n{\F2 fj Fl\)n-l^m,n 
+n+l{\F^fiFl\)Jm,n+l + n-l^ fiFl\)Jm,n+l] + • • • , (17) 

where the ellipsis reflects terms which are vacuum on-diagonal. Use is made of the fact 
that both and the projectors p| = fjf^ andp^ = /j// commute with the Fn{{Snz}) 
and each other, and that e.g., ujjfi = N-^'^ Y.kk'k" e*'' *'"e*'='-^"e'*^" *'"ti-/l /-,. The 
vacuum ofi^-diagonal part corresponds to contracting out a factor of S^/N. This occurs 
in two ways for each term to give (17). The two vacuum diagonal terms have a different 
number of particles in the final state, corresponding to (spin) hole and and particle 
excitations. The result can be written symbolically as 

Ui = lifiS^ + flS- + S+f, + S-fh + ..., (18) 

with the understanding that for, e.g., in 5+/^ the 5+ acts only on the |)o vacuum 
to the left. The ellipsis represents terms which are necessarily diagonal in the bare 
vacuum state. Since m„ is unitary Uj^\S) is a normalised state vector and it is to be 
noted that the four vacuum off'-diagonal terms displayed above exhausted exactly one 
half of this normalisation. The on-diagonal ellipsis, evidently, has equal weight. 

That, when it acts on the physical vacuum, the flux tube operator can create 
excitations (with a cancelling change in the total Sz) is also evident from a comparison 
between Szn + 1/2 = and 5+ = /^m„. The effect of these two operators acting 
upon the singlet physical vacuum must be essentially the same, i.e., the effect of 
w„ is the same as /„ but without a change in Sz- Multiplying (18) by /+ gives, 
Sn = flfniS^/N) + {S+/N)fjX. (There is no commutator with 5+ in the last term 
because of the understanding that this acts on the vacuum to the left.) 

Upon examining the ensemble of physical operators it should be observed that 
these can be consistently interpreted in terms of /-particles which are S = \/2 spinous 
with Sz = +1/2 while the flux-tube-particles are the equivalent with Sz = —1/2. 

For the t — J model, in the absence of doping, the spin sector kinetic energy is 
generated by the vacuum diagonal part of the transverse exchange, i.e., 

+ (19) 

For a near uniform state this must have an expansion +4 fJ uaiI f, ~ ci^fjfi + 
fifj S CL^n^'^m whcrc thc Sum arc deviations from the imiform state. Since, by defini- 
tion a} = {u^u}j) does not depend upon the uniform state involved, it can be evaluated, 
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with a suitable gauge, for a state veetor which is an equal weight of the two Neel states. 
Both states have a uniform field with one half a flux quantum per plaquette and in a 
(non- JW) gauge which reflects the symmetry of the lattice this implies that the change 
of phase introduced by the UjuJ — > e^"/"*. There are no direct matrix elements be- 
tween these two state and so is givc^n by a simple average of these two phases, i.e., 
the leading term site diagonal part of (19) is, 

+ ^cosJ/t/^.^:^/t/,, (20) 

with J' = J I\f2. Thus, e.g. for a two dimensions the bare energy of an excitation of 
momentum fc is eg = 2J'7g. 

The next term in the expansion of (19), written as -\-^SfS'^ , is obtained using the 
vacuum off-diagonal part of 5+ deduced above. There are four equivalent contributions 

to give the result, 

+Jflfifjfj + H.c. +Jflf]{fjh) + H.c. = ±Ao///j + H.c, (21) 

admits a pair amplitude. With suitable signs, the Fourier transform corresponds to 
rf-wave pairing, i.e., a term, 

^EW.^ + ^-' (22) 
fe 

where A^: = Ao(cos/;;x — cos ky). Clearly, the original exchange term does not change 
the value of Sz and the pairing in (22) is to be interpreted in terms of k, Sz = 1/2 
and —k, Sz = —1/2 spinon pairs. 

The final version (15) of the JW formalism would appear very similar to the usual 
auxiliary (slave boson) particle scheme. The down spin fermion has been replaced by 
a flux tube operator However this parallel is very misleading. In terms of up and 

down fermions the transverse exchange is: +^ J2<ij>{f}ifiifljfij +H.c.). Invariably 
in the "slave boson" approach this interaction is subjected to a mean field factorisation: 

-Ei(/iVii)/TVTi + ^^-'=-' (23) 

<ij> 

which changes the sign and reduces the magnitude of the matrix elements for the 
hopping of the spin particles. This latter philosophy is very attractive since it separates 
the up and down degrees of freedom and leads to an identification of the /j^j as spinons, 
i.e., spin 1/2 particles. The cost is that the constraint is all but sacrificed. The point 
of the present formalism is precisely to maintain this constraint exactly and yet to 
show that the excitations arc indeed spinons. 

The present formalism has its own naturel approximations. If t :» J, the motion 
of the holes is certainly rapid compared to the dynamics of the lowest lying spin levels. 
If the concentration of holes is sufficient there will occur a motional averaging of the 
spin wave function. The co-movement of the spins, as the holes hop, will cause the 
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spinon-spinon interactions to be averaged and suppress the associated corrcilations 
(and ordering). Given their rapid motion, the holons can be added by perturbation 
theory. On the charge time scale the spins can be considered as frozen. Consider a 
given up spin at site n, it sees itself in a virtual crystal in which there is a certain 
density of holes at near neighbour sites. This up spin then admixes with the holes on 
the neighbour sites and this implies a probability ~ a;f/|/z| ^ x of observing a hole at 
n. Whatever the details, it is implied that the state which was fl\) now becomes, 

{afl + pbiui)\)^Pj), (24) 

where |ap + = 1 with the concentration of holes x = and where the flux 
tube is needed in order that have fermion commutations rules. A Fermi sea is 
now constructed from the Fourier transform of the The combination found in 
implies a pairing between the 6jj holons and /„. This can be associated with Kondo 
physics since with a small (3 the excitations are relatively heavy. Likewise, the state 
for a down spin site, 

ui\)^{aui + l3biui)\). (25) 

With Sz = the presence of a //^ particle implies a probability of 1 — a; of finding 
an up spin at site n, and there are N/2 such particles independent of x. In total, 
there are {N/2){1 + x) degrees of freedom for such a system and the present "Kondo" 
approximation implies that there are only N/2 low energy energy excitations and that 
these are predominantly "spinon-like" with a small admixture of holon character. For 
low enough energies, the remaining holes are reflected, (25), as a renormalisation of 
the vacuum. The site independence of /? implies a holon condensate with q ^ 0. The 
admixture of holon character into the /-particles implies that the term, e.g., —tblu^u'jbj 
contributes to the hopping matrix elements for these renormalisation particles. As a 
result J' — > Je = (1 — x)J' — 2xt. With this approximation, the effective Hamiltonian 
is, 

k k 

All interactions between spinous and holons which appeared in (15) have been dropped 
since these will be average towards constants by the rapid hole motion. The physical 
conduction electron propagators (Itc|j|(t)c^|(0)) = c{Tt f ^{t) f}^{Q)) , i.e., are c, the 
condensed faction of holons, times the spinon propagators. 

Of considerable importance is the fact that the kinetic energy is proportional 
to Je = (J' — 2xt). Nagaoka's theorem [3] implies that a small hole doping favours 
ferromagnetism. The ferromagnetic, i.e., negative sign of the 2xt reflects this fact. 
(This is the itinerant hole version of ring exchange.) Clearly the anti-ferromagnetic 
J'-term must have the opposite sign. As a consequence, in absence of pairing^ this 
Hamiltonian would exhibit a quantum critical point (QCP) at Xc = {J' /2t) at which 
the spinon band collapses. Clearly the specific heat and magnetic susceptibility diverge 
at such a point. Since superconductivity is favoured by such a large density of states, 
will be a maximum near Xc and superconductivity will "hide" this QCP. This might be 
described as a pseudo-QCP with the real QCP being shifted to smaller concentrations 
at which the holons condensate and superconductivity first appears. 
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